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A Search for Star-Disk Interaction Among the Strongest X-ray Flaring Stars in 

the Orion Nebula Cluster 

Alicia N. Aarnio^, Keivan G. Stassun^^^, and Sean P. Matt'^ 
ABSTRACT 

The Chandra Orion Ultradeep Project observed hundreds of young, low-mass stars under- 
Koing highly energ etic X-ray flare events. The 32 most powerful cases have been modeled by 



Favata et all (J2005l) with the result that the magnetic structures responsible for these flares can 

be many stellar radii in extent. In this paper, we model the observed spectral energy distribu- 

^N , tions of these 32 stars in order to determine, in detail for each star, whether there is circumstellar 

^^ .... . . ... 

disk material situated in sufficient proximity to the stellar surface for interaction with the large 

magnetic loops inferred from the observed X-ray flares. Our spectral energy distributions span 

the wavelength range 0.3-8 /im (plus 24 /xm for some stars), allowing us to constrain the presence 

of dusty circumstellar material out to > 10 AU from the stellar surface in most cases. For 24 

i^ ■ of the 32 stars in our sample the available data are sufficient to constrain the location of the 

I , inner edge of the dusty disks. Six of these (25%) have spectral energy distributions consistent 

O I with inner disks within reach of the observed magnetic loops. Another four stars may have gas 

"^ ' disks interior to the dust disk and extending within reach of the magnetic loops, but we cannot 

d . confirm this with the available data. The remaining 14 stars (58%) appear to have no significant 

disk material within reach of the large flaring loops. Thus, up to ~ 40% of the sample stars 

exhibit energetic x-ray flares that possibly arise from a magnetic star-disk interaction, and the 

remainder are evidently associated with extremely large, free-standing magnetic loops anchored 



only to the stellar surface. 



^^ I Subject headings: stars: circumstellar matter, stars: flare, stars: pre-main-sequence 

o 
o. 

1. INTRODUCTION 

The recent large x-ray surveys of the Orion and Taurus star-for ming regions performed by Chandra 
rS ' and XMM (i.e., COUP, XEST: lGetman et al.ll2005aHAudard et al.ll2007 ) provide an unparalleled opportunity 



jrt ' to study the magnetic activity of young, low-mass stars. These deep observations spanning long temporal 

baselines (e.g., the coup x-ray light curves span 13 days with near-continuous time coverage) reveal that 
low-mass pre-main-sequence (pms) stars possess x-ray luminosities 3-4 magnitudes greater than that of the 
present-day Sun and exhibit extremely energetic flaring events with high frequency. 

Detailed analyses of these flares reveal that they are simi lar to solar fl a res, b ut are orders of magni- 



tude more energetic and larger in physical size. In particular, iFavata et al.l (|20Q5h subj ected the 32 most 



energetic flares observed by COUP to analysis via a standard uniform cooling loop model (jReale et al.lll997 : 
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Svlwester et al.lll993l : iPriest fc Forbesll2002t iFavata fc Micelall2003l ). with which they derived the properties 
of the magnetic coronal loops that participate in the flare events. They found that these magnetic loops 
were extremely large — extending tens of stellar radii in some cases — much larger than ever observed on older 
stars. Such large-scale flares could have important ramifications for a number of issues, such as the shedding 
of stellar angular momentum and mass, powering of outflows, and ionization/dissipation of circumstellar 
disks. 

Magnetic loops with sizes on the order of ~ 10 stellar radii have long been postulated as part of 
magnetospheric accretion scenarios. In this paradigm, a large-scale stellar magnetic field th r eads the inner 
edge of a circumstellar disk, channeling accretion fr om disk to star (e.g., 



Shu et al.lll994ilHartmannlll994tlHavashi et al.lll996[ ). Indeed. iFavata et al. 



Camenzindlll99d lKoenigill991: 



( 20051 ) speculated that the large 



magnetic loops observed in the COUP sample may be facilitating this type of magnetic star-disk interaction, 
in part because they argued that such large loops would likely be unable to remain stable if anc hored only 
to the stellar surface. However, as most of the 32 COUP sources studied by IFavata et al.l (|2005l ) lacked at 
that time sufficient photometric data with which to characterize the optical-infrared (ir) spectral energy 
distributions, they could not confirm the presence of inner disks to which the observed magnetic loops might 
link. 

Thus, there is still an outstanding question as to whether stellar coronal activity in these stars alone 
can drive such energetic flare events, or whether the energy (or at least the trigger) derives from a star-disk 
interaction. While the l atter requires magnetic loops large enough to reach the inner edge of the disk, several 



theoretical studies (e.... lu.denskv et al.l Ei lMatt X. Pudrit^El and references therein) have shown that 



the presence of a disk truncates the stellar magnetosphere so that closed magnetic loops extend not much 
further than the inner edge of the disk. Thus, if these flares are powered by the star-disk interaction, one 
would expect the size of the fiaring loops to approximately coincide with the location of the disk inner edge. 
On the other hand, if the energetic flares are purely a stellar phenomenon, the largest loop sizes may only 
be exhibited by stars that lack disk material close to the star. To address this question, it will therefore be 
useful to be able to determine the proximity of disk material to the furthest extent of the flaring magnetic 
loops exhibited by the sample analyzed in this work. 



Near- and mid-iR colors can be used as a crude tracer of close-in circumstellar material. Getmanetal 



(J2008bl ) have used Spitzer IR colors to distinguish Class ii and Class in objects (i.e., stars with dusty disks 
and naked T Tauri stars, respectively) among 161 flaring COUP stars. Interestingly, they found evidence 
that whereas the largest flaring loops tended to be associated with Class III objects, the Class ii sources in 
their sample were more likely to possess relatively small magnetic l oops. The IR co lors alone do not provide 
a quantitative measure of the location of the inner disk edge, but ICetman et al.l (12008b) suggest that the 
magnetic loops may be confined by the disk to be within the disk co-rotation radius (the radius at which disk 
material, if present, orbits the star with angular velocity equal to the star's angular velocity). This suggestion 
is important, as the disk co-rotation radius is the point spec ifically at which some magnetospheric accretion 
theories predict magnetic star-disk interaction to occur (e.g.. lOstriker &: Shulll995[) . It is desirable, therefore, 
to establish the relationship between circumstellar disks and the very large magnetic loops observed by COUP 
more quantitatively than IR colors alone permit. 

In this paper, we present detailed spectral energy distributions (SEDs) for each of the 32 most powerful 
x-ray flaring COUP sources at wavelengths 0.34-8 /xm, plus upper limits at 24 ^m ([J2]). In SjSjwe present 
near-infrared color excesses for the sample as a basic tracer of close-in circumstellar disks. Next we compare 
in detail the full observed SEDs against synthetic SEDs of low-mass pms stars with disks (§||4]'[5]) in order 
to (a) ascertain whether dusty disks are present around these stars, and (b) if so, determine quantitatively 



whether the inner edges of those disks are su fhciently close to the stellar surface to interact with the large 
flaring loops observed by iFavata et al.l (|2005l ). The results (SJH]) indicate that approximately more than half 
of the sample stars lack significant disk m aterial within r each o f their flaring magnetic loops; evidently the 
extremely large flaring loops observed by iFavata et alj (J2005l) are in most cases free-standing structures 
anchored only to the stellar surface. In ij7]we discuss some implications of this finding. 



DATA 



2.1. Study Sample, Loop Heights, and Stellar Data 



The 32 stars for our study constitute a unique subset of the COUP (|Getman et al.ll2005bl) observation, 
identified by IFavata et al.l ( 20051 ) as exhibiting the brightest r-^ 1% of all flares observed by COUP. These 32 
flares had sufficient photon statistics with which a uniform cooling loop (ucl) analysis could be performed. 

The UCL model is based on observations of solar flares. The occurrence of reconnection events on the Sun 
has been used to benchmark relationships between x-ray flare decay slopes and the magnetic fleld structure 
conflning the emitting plasma. After a magnetic reconnection event occurs, heated plasma evaporates from 
the c hromosphere into the confining loop. The material then emits soft x-rays as it cools ( Priest fc Forbes! 
20021 ). and the x-ray light curve's decay tiru e as well as its slope in density-temperature space is related to 
the magnetic loop length (JReale et al.lll997l ). This method was developed using hydrodynamic simulations 
which were calibrated against spatially resolved imaging observations of solar flaring loops. 

Favata et al.l (|2005l) applied the ucl analysis to their sample of 32 stars and t hus derived the lengths 



of the magnetic loops confining the flare events observed by COUP. A reanalysis by iGetman et al.l (|2008al ) 
includes these 32 objects. The derived loop lengths in both studies are consistent within uncertainties, so 
we adopt the former for consistency throughout. For simplicity, we estimate the loop height from the stellar 
surface as the loop length divided by 2. This is an upper limit; for example, in a circular or semi-circular 
loop geometry, the actual loop height would be the loop length divided by n. These loop heights (half loop 
lengths) and their uncertainties are summarized in Table |TJ Uncertainties in the loop lengths are due to 
uncertainty in the measurement of the flare's peak temperatur e, decay time , and d ecay slope. For a detailed 
discussion of the quoted uncertainties in the loop lengths, see IFavata et al.l ( 20051 cf. their §3.2). 



To narrow the range of acceptable best-fit spectral energy distributions {^, we require basic stel- 
lar paramet e rs inc luding effective temperatures (Tog), and radii (iistar)- These are taken primarily from 
Hillenbrandl (|l997l) and are summarized in Table [T] In cases where stellar parameters were not available 
from the literature, we adopt the temperature s and radii of the best- fit SED model (see S|4|). Table [T] also con- 
tains Ca 11 equivalent widths (as measured by iHillenbrandll 1 9971 ) and our newly reported AKs and A{U — V) 
excess measurements (see Sj3]) . We use AKs excess as a supplemental indicator of close-in, hot disk material, 
and the Ca II equivalent widths in combination with A(U — V) excess to indicate ongoing accretion onto the 
stellar surface. 



2.2. Photometric Data 



Fluxes for each of the 32 stars in our study sample were assembled over the wavelength range 0.34 /im 
(u band) to 24 /im {Spitzer Multiband Imaging Photometer, mips). With these data, we probe the stellar 
photosphere and circumstellar dust content. These measurements are summarized in Tables [5] and [31 



Optical fluxes were taken from the ground-based observations of the Orion Nebula Cluster (ONC) by 
Da Rio et alj ( 20091 ) in the UBVIc passbands (0.36, 0.44, 0.55, and 0.83 /im, respectively), obtained with the 
ESO Wide Field Imager (wFi). We supplemented thes e with fluxes from Hubble Space Telescope Advanced 
Camera for Surveys (acs) data (JRobberto et al.l 120051) . providing broadband flu xes at 0.43, . 54, .77, and 
0.91 /im, as well as V and Ic magnitudes from the ground-based observations of lHillenbrandl (J1997I ). 

The 2 Micron All Sky Survey (2MASS: ISkrutskie et al.ll2006l) provides near-infrared JHKs magnitudes. 
Of critical importance to our analysis, infrared photometry from the Spitzer Infrared Array Camera (iRAC) 
and MIPS instruments provide the clearest probes of warm circumstellar dust in the 3.6, 4.5, 5.8 and 8 /im 
bandpasses (irac) and at 24 /im (mips). We measured these fluxes using pipeline-processed, archival data, 
and found our values to agree within a few mJy of the unpublished measurements of the Spitzer GTO team 
(S. T. Megeath, private communication). The mips fluxes were measured by us from the Spitzer archive 
using the pipeline reduced 24 /im images. Unfortunately, the mips image of the ONC is saturated over most 
of the region of interest, and we were thus unable to recover more than a few upper limits (see Table [3]). 

Where magnitudes were originally reported, these have been converted to fluxes using published zero 
points for each instrument (see flnal rows in Table [5] and [3]) . In addition, we have in general adopted larger 
uncertainties on the fluxes than the formal measurem ent errors, in order to account for typical v ariability lev- 
els in the optical and near-infrared of ~ 0.1 mag (e.g.. lllerbst et al.lll994HCarpenter et al.ll200ll ). Speciflcally, 



we adopt an uncertainty of at least 10% on the fluxes, unless the formal measurement error is larger. 



PRELIMINARY DISK DIAGNOSTICS: COLOR EXCESS AND ACCRETION 

INDICATORS 



Traditionally, the presence of warm circumstellar dust around low-mass pms stars has been traced using 
near-iR "color excesses," such as A{H — K), deflned as the differe nce between the observed (de-reddened 
color and the color expected frorn a ba re stellar photosphere (e.g.. lStrom et al.lll989t iLada fc Adamslll99! 



Edwards et al.lll993l : iMeyer et al.lll997l ). The use of a single color excess of course does not permit a detailed. 



quantitative determination of disk structure (such as the size of the inner truncation radius, which is our 
primary interest here) because a given color excess depends in complex ways upon multiple disk and stellar 
parameters (see also Sj4lbelow). Moreover, near-iR colors may cause the observer to miss the presence of some 
disks, particularly those with large inner holes or around very cool stars whose photospheres peak in the 
near-iR. For example, the overall disk frequency in the Orion Nebula Cluster has been estimated at ~65% on 
the basis of excess emission in the K band (Hillenbrand et al. 1998), but increases to '^85% simply by adding 
an L-band measurement (Lada et al. 2000). In other words, the addition of mid-iR measurements can be 
very important for the detection of disks ( for a discussion of dis k-detection efficiency using near-iR colors, see 
Hillenbrand et al.lll998t iLada et al.ll2000t lErcolano et al.ll2009l ). This is particularly relevant to our sample, 
in which many of the stars have cool photospheric temperatures, and where even disks with relatively large 
inner holes could be within reach of the observed very large magnetic flaring loops (see Table [T]) . 

Still, color excesses have the advantage of being easy to collect and analyze for large numbers of stars — 
especially prior to the advent of the wide and deep longer-wavelength surveys made possible by Spitzer — 
and of providing a relat ively straightforward "yes/no" criterion for the presence of a disk. For example, 
Hillenbrand et al.l (jl998l ) used the A(/ — K) color excess to conduct a census of disks among ~ 1000 low- 
mass stars in the ONC. In that approach, the observed V — I color was used to measure the extinction. Ay, 
which was used in turn to deredden the observed I — K color. Any excess A{I — K) was then attributed to the 



presence of a disk. This approach has the advantage of requiring for each star only three flux measurements 
{VIK) and a spectral type (with which to establish the expected photospheric colors). It assumes that (a) 
the observed /-band flux is purely photospheric in origin, and (b) the observed V-band flux is only affected 
by reddening (i.e., does not include any "blue excess" due to veiling emission from accretion). 

Usin g our compiled fluxes in T ables [2] and [3] for our study sample, we have calculated A.Ks in a manner 
similar to iHillenbrand et al.l (|l998l ). However, rather than use only the observed V — I color to determine Ay, 
and rather than normalize the stellar flux to the observed /-band flux, we have performed a two-parameter fit 
to each star's observed SED. To isolate the stellar flux from the disk and/or accretion flux, we use fluxes which 
appear to be photospheric in ori gin only, excluding th e bluest wavelength fluxes which could be affected by 
accretion flux or scattered light (jWhitnev et al.ll2003al ) as well as the reddest wavelength fluxes which could 
contain flux from a disk. In general we used the fluxes from 0.5 to I.O /im (total of 8 flux measurements, see 
Tables [2] and [3]) for this fitting, which should be a substantial improvement over the two- band (V and Ig) 



appro ach described above. To the observed fiuxes we fit a NextGen model atmosphere (jHauschildt et al 



19991 ) at the spectroscopically determined Tcs from the literature (Table [T]). The two free parameters of the 
fit are the Ay and the o verall normalization of the stellar flux. For comparison, previously published Ay 
values (|Hillenbrandlll997l ) are listed in Table [TJ and our newly determined Ay values are also reported along 
with their corresponding errors (99% confidence limits) as determined from the two-parameter SED fit. For 
10 stars in our study sample. Ay values are reported here for the first time. In ma ny of these cases we find 
large Ay values {Ay > 10). These stars were likely absent from the optical study of lHillenbrandl (|l997[ ) due 
to the high reddening/extinction. 

For about half of the stars in our sample, our newly measured Ay values agree within 99% confidence 
with those previously reported. For the remaining stars, the Ay values differ significantly; the reason for 
this difference is illustrated in Fig. [1] for the case of star COUP 262. The red model SED re presents the 
previo usly reported fit of a Tcs = 4395 K photosphere to just th e V and Ig band fl uxes from 



Hillenbrand 

()l997l ). The extinction that results is Ay = 3.77 (see Table fTJ and lHillenbrandlll997l ). an d as a consequence 
the / C-band flux appears to be highly in excess of the photosphere with A(/ — K) = 2.21 (IHillenbrand et al 



19981 ). However it is evident from visual inspection of the complete set of observed fluxes that this model flt 
is a poor representation of the additional measurements included here. As shown in the flgure, our new flt 
to the entire set of fluxes gives Ay = 7.91^° 53 and AKs = 0.14 ± 0.11. For this particular case, evidently 
the previously reported y-band flux was anomalously high by ^ 5a (perhaps due to the ubiquitous optical 
variabflity of PMS stars), and so fltting for Ay to just the V and Ig fluxes resulted in a distorted model SED. 
The difference between the two SED fits is important in the context of our study: e.g., the previously reported 
value of A(/ — A') for COUP 262 implies a massive, warm circumstellar disk close to the star, whereas our 
newly determined best-fit SED model is in fact consistent with no close-in disk. We revisit the SED fit of 
COUP 262 in the context of the entire study sample below (§§|4]and[6]). 



Hillenbrandl (|l997l ) adopted photospheric colors of main-sequence dwarfs (i.e., for dwarfs, log 5 ^ 4.5) 
in the calculation of near-iR excesses. However, low-mass PMS stars at the young age of the ONC (^ 1 Myr) 
are expected to have somewhat lower log 5 valuesdue to their large radii. Thus we have also considered the 
extent to which the assumed log 5 affects the predicted stellar colors and thus the inferred near-iR excesses. 
In Fig. [5J we compare the SEDs of NextGen stellar atmosphere models as a function of log 5 for several 
representative Tag appropriate for our study sample (all models shown are solar metallicity). We see that 
the choice of logg is not important for stars warmer than T^e ^ 4200 K. However, we find that for cooler 
objects the predicted I — K colors become increasingly redder with decreasing log 5, which could lead to 
over-estimated near-iR excesses for the coolest stars. Therefore in this study we have opted to use model 



atmospheres with logg appropriate to each object {logg calculated from previously reported stellar masses 
and radii; see Table [T|). 

Newly determine d AKs values following the procedure described above are reported in Table [J Several 



previ ous studies (e.g. JStassun et alJll999HRebullll200lUHerbst et alJbQOa iLamm et al.ll2004 iMakidon et al 



20041 ). adopted a threshold value of Ai^s > 0.3 (i.e., ~ Scr excess given typical aK = 0.1 mag; see Sec. 
for identifying stars with close-in circumstellar disks. By this criterion alone, five of the stars in our sample 
(COUP 141, 223, 1246, 1343, 1608) show large near-iR excesses indicative of the presence of warm circum- 
stellar dust 1^. The remaining stars in our sample show very weak or no evidence for near-iR excess emission 
(i.e., AKs < 0.3). 

In the context of the principal aims of the present study — where we seek to determine whether the 
large magnetic loops observed in the sample stars are linked to circumstellar disks — one might anticipate 
that AKs could be used as a quantitative tracer of circumstellar dust located within reach of the observed 
magnetic loops. For example, if the presence of substantial AKg excess correlates with the location of the 
disk truncation radius (i?trunc: the distance from the inner edge of the disk to the star), then we might simply 
take stars with AKa, > 0.3 as those whose i nner di sks are likely to be magnetically linked to the star. This 



is similar to the approach of iGetman et al.l (|2008bl ). 



The combination of factors discussed above, however — relatively cool stellar photospheres which peak 
in the near-iR, magnetic flaring loops that are large enough to interact with disks at relatively large heights 
above the stellar surface, etc. — makes near-iR excess an inefficient tracer of the types of disks we seek to 
characterize. Table |4] categorizes the sample stars according to whether or not they display significant near-iR 
excess emission (i.e., AKg > 0.3 vs. AKs < 0.3) and whether the inner-disk truncation radius is larger or 
smaller than the dust destruction radius (i.e., i?trunc > -Rdust vs. i?trunc ^ ^dust, where -Rdust is the distance 
from the star within which dust is warm enough to sublimate) . Here -Rtrunc is determined from our detailed 
SED model fitting as described below (§4]). The off-diagonal elements of Table |4] represent cases contradicting 
the assumption that AKs excess correctly and quantitatively predicts the location of the inner-disk edge. 
For only one of the sample stars (COUP 1246) do we find a relatively large inner-disk hole (i?trunc > -Rdust) 
but a large AKs > 0.3. This one case appears unusual because the best-known explanation for a strong 
near-iR excess is the presence of warm dust close to the star. Upon closer inspection of this case (COUP 
1246), we found that the fitted Ay was underestimated because of the presence of a moderate blue excess, 
likely due to chromospheric activity as suggested by the observed filled- in Ca II emission (Table [l|) , and 
which is not included in the photosphere model. If we manually adjust the Ay value upward by '~ 3cr from 
the fit value of Ay = 1.52to:6o (see Table[l|) to Ay = 2.22, the value of AKs becomes 0.17. Indeed, our final 
best fit SED model (see ^ has an Ay of 3.04, and thus an even lower AKs- Therefore, for the following 
analysis and discussion, we assume that the calculated excess AKs in COUP 1246 is not significant (i.e., 
AKs is consistent with being less than 0.3). 

Several stars in our sample exhibit a meager AKs < 0.3, and by this criterion alone would be classified 
as lacking close-in, hot dust. However, for eight of these stars (in the lower left quadrant of Table 4), our SED 
fitting found these to in fact have disks that reach relatively close to the star (i?trunc ^ -Rdust), and in some 
cases (shown below, ^, the disk reaches sufficiently close to the star to interact with the observed large 
flaring loops. In summary. Table 4 indicates that while AKs ^ 0.3 appears to be an accurate indicator of 
dusty material close to the star, the lack of significant AKs does not rule out the presence of dusty material 



^As described below, the large AKg in COUP 1246 is likely the result of an underestimated Ay, and we consider its near-lR 
excess to be not significant for the remainder of our analyses. 



close to the star. Thus, for the purposes of the present study where we seek to estabUsh more quantitatively 
the location of the inner-disk edge in relation to the observed large magnetic flaring loops, we cannot rely 
solely on traditional near-iR excesses. 

In addition to AKs measurements, from our new SED fits we calculate blueward color excesses, A{U — V) 
which can be used as a tracer of accretion at the stellar surface ( "hot spots" ) and/or chromospheric activity 
and can help in our interpretation of some SEDs (see ^ these are reported in Table [T]). Table [S] categorizes 
the sample stars in a manner similar to TablelH but now using A{U — V). We i dentify objects with A{U — V) 
< —0.3 as those likely possessing hot accretion s pots on their surfaces (se e Rebull et al.l 120001 ) and thus 
likely to be undergoing active accretion (but see iFindeisen fc Hillenbrand! l2010l for a discussion of other 
phenomena that may cause blue excesses in pms stars). Five of the sample stars show evidence for active 
accretion, and all but one of these have i?trunc ^ ^dust as expected for a disk that extends clo se enough to 
the st ar for accretion to occur. Furthermore for two of these stars the Ca II measurements of iHillenbrandl 
(jl997t see Table[T]) also indicate active accretion. Only one star (COUP 1568) fails to show A{U — V) excess 
despite possessing a close-in disk edge. Thus, while A{U — V) cannot provide a quantitative measure of the 
location of the inner disk for non-accretors [A{U — V) > — 0.3] , it is as expected a relatively reliable indicator 
of ^trunc ^ ^dust for active accretors [A{U — V)< —0.3]. 

In the analysis that follows, we use detailed SED fits over the full range of available photometric data 
(Tables [2] and [3]). Where applicable, we use the AKs and A{U — V) excesses and Ca II equivalent widths 
in Table [1] as secondary information to aid our classifications in order to characterize in detail the presence 
and structure of circumstellar disks in our sample. 



SYNTHETIC SPECTRAL ENERGY DISTRIBUTION MODELS 



To compare the observed SEDs of our sample with the SEDs expected from young stars wi th disks within 



reach o f the observed flaring loops, we employed the Monte Carlo radiative transfer code of I Whitney et al 



(|2003a| bl). ttsre, to generate synthetic SEDs. Our aim is to more quantitatively constrain the structure of 
any circumstellar material around each of the 32 stars in our sample so that we may determine , in detail for 



each s tar, whether there is in fact disk material within reach of the magnetic loops observed bv lFavata et al 
(J2005l ). Thus for each of the stars in our sample, we wish to determine the range of disk parameters — the 
most important of these being the location of the inner edge of the disk — that are able to reproduce the 
observed SEDs {^. 

The TTSRE code models randomly emitted photons from the central illuminating source and follows 
the photons as they interact with (i.e., are absorbed or scattered by) any circumstellar material. The 
circumstellar material is modeled as an optically thick dust disk extending from an inner truncation radius, 
-Rtrunc, to an outer radius of typically a few hundred AU. The disk in general may be "flared" such that its 
scale-height increases with increasing distance from the star, or it may be flat. Surrounding the star and 
disk may be a spherically distributed infalling envelope with bipolar cavities; such an envelope is generally 
required for reproducing the scattered-lig ht properties of embedded objects (generally seen as moderate 
excesses in the blue; e.g.. IStark et al.ll2006f ). The code also self-consistently solves for thermal equilibrium 
in the disk as absorbed photons heat the disk and are re-radi ated. Subhmation of du st is also included (for 
details of the dust properties used by the code, see Table 3 of I Whitney et al.l 12003 ol. The co de models the 
central illuminating source using the NextGen atmosphere models of lHauschildt et al.l (jl999l ) . We adopted 
the solar- metallicity atmosphere models, with logg and Toff chosen according to each star's observationally 



determined Afgtar, Rstm-, and Tcfi (Table [T]). 



Because of the very large number of permutations on the possible star/disk/envelope parameters in- 
cluded in the ttsre model (i.e., disk mass, disk inner and outer radius, disk flaring profile, disk accretion 
rate, disk inclination angle, etc.), there is in general not a simple one-to-one correspondence between a given 
observed SED and, say, i?trunc- Thus, to fully explore the range of disk parameters that could possibly repro- 
duce the observ e d SED s of our sample, we made use of the very large grid of ttsre models constructed by 
Robitaille et al.l ( 2006[l . The grid in cludes some 200,000 models representing 14 star/disk/envelope parame- 
ters (see Table 1 of lRobitaille et al.ll2006) that were independently varied to encompass virtually all possible 
combinations for young stellar objects with masses 0.1-50 M0 in the Class O-iii stages of evolution. The 
parameter space for this grid was specifically set to be very finely sampled for T Tauri stars (i.e., TcS <5200 
K); at higher temperatures, >;5200 K, the grid is more sparsely sampled. Whereas at TcS <5200 K the grid 
is sampled in i?trunc by ^5%, at Toff > 5200 K it is sampled much more sparsely at ~50%. This issue affected 
only the hottest object in our sample, COUP 597 (see ^^. Additionally, by construction, the models in this 
grid are set to only include non-zero accretion rates and to always include emission from a hot accretion spot 
on the star. As a result, in some cases blueward excess due to the hot accretion spots is seen in the best- fit 
model SEDs, which we disregard when we lack U or B band fluxes to constrain the blue side of the SED. 

With the added free parameter of extinction. Ay , we searched the grid via x^ minimization for all syn- 
thetic SEDS that fit the observed SED of a given star within the 99% confidence level (that is, we re iected those 
models that yielded a Ax^ goodness-of-fit likelihood of 1% or less relative to the best- fit model: iPress et al. 
19951 ). For stars with spectral- type determinations from the literature (see Tabled]), we furthermore require 
the model fits to have stellar T^g within 500 K of the literature value except for a few cases where we found 
it necessary to relax the Tcfi constraint in order to achieve an acceptable SED fit; these exceptions are noted 
when we discuss each object individually below. 

In addition to temperature, we also filter the best-fit models by disk mass. In fitting the SEDs of the 
sample stars, particularly in cases with little or no IR excess emission in the observed SED, we found that a 
number of the best- fit models none theless had disks with s mall Rtrunc, but only if the disk also had a very low 
mass. The extensive model grid of lRobitaille et al.l (|2Q06r ) allows for disks with masses as low as 10"^" M©. 



Such low disk masses, however, may be well below what is physically realistic, and certainly below what 
is observable, for young T Tauri stars. Recent detailed studies of p ms stars with so-called "transitional" 
and "pre-transitional" circumstellar disks fe.g. . lEsoaillat et al.ll2007l) — disks that are undergoing the rapid 
disk-clearing process from the inside out (e.g., iBarsonv et al.l 12005) — show that even at this late stage the 
circumstellar disks are in fact quite massive. For example, lEspaillat et al.l (|2007l ) derive Mdisk ~ 10~^ Mq 
for the pre-transitional disk of LkCa 15 with i?trunc ~ 45 AU, and Mdisk ~ 10~^ M© for the slightly more 
evolved disk of UX T au A, with -Rtrunc ^ 60 AU. A more extreme case is that of CoKu/Tau 4, for which 
D'Alessio et al.l (|2005r ) find an extremely low Afdisk ~ 10^^ M0. In what follows, we will thus restrict our 
analysis to include only model SEDs with Mdisk > 10^^ M0 as more accurately representing the empirical 
disks of young, low-mass stars. For illustrative purposes, however, we display all models with disk masses 
greater than lO""* Mq. 

In addition, for the purpose of interpreting the resulting best-fit model SEDs, we found it useful for each 
star to generate an additional ttsre synthetic SED as a fiducial reference model. For the cases where there is 
apparent IR excess in the data indicating the presence of a disk, we generate a fiducial model identical to the 
best-fit SED model, except that we set -Rtrunc equal to the magnetic loop height, -Rioop (TableU]). In the cases 
for which the data show no IR evidence for dusty disks, the fiuxes of the best fit SED model are essentially 
arbitrary beyond the longest wavelength data point. For simplicity in these cases our fiducial model is a 



simple star+disk SED, adopting stellar properties from the literature. The modeled structure is that of an 
optically thick, geometrically thin, slightly flared disk with no envelope and no accretion. The disk mass in 
these cases is set to 0.01 Mq, and its inner truncation radius is again set equal to the magnetic loop height. 
Thus, in all cases, the fiducial model allows a direct, visual comparison of the observed and best-fit SEDs 
against that expected if the inner disk is within reach of the magnetic flaring loop. It is important to note 
that in cases where the magnetic loop height is within the dust destruction radius, the disk is truncated at 
dust destruction by default (i.e., the ttsre models require i?trunc > ^dust)- 

We present in this work the SEDs of the 32 stars in our study sample. Flux measurements and upper 
limits (Tables [2H3]) are represented by diamond symbols with error bars (which are in most cases smaller 
than the symbols) and blue triangles, respectively. The SED of the underlying stellar photosphere (NextGen 
atmosphere) is shown for co mparison as a d a shed line. Superposed on the observed SEDs, the best-fitting 
SED models from the grid of iRobitaille et al.l ( 20061 ) discussed above are shown as dash-dotted curves. The 



fiducial TTSRE model that we calculated is shown as a solid, blue curve. Figures [SHni are included in print 
for illustrative purposes below (see ^, with the remaining figures accessible electronically. 



5. INTERPRETING THE SPECTRAL ENERGY DISTRIBUTIONS 

With observed and model SEDs in hand for all 32 stars in our study sample, we can now attempt to 
answer the central question of this paper: Are the large flaring loops observed on these stars likely due to 
a magnetic star-disk interaction, or do they represent primarily stellar phenomena unrelated to disks? To 
answer this question, in this section we discuss the specific criteria by which we determine, from examination 
of each star's SED, the likelihood that it possesses a disk whose inner edge is within reach of the observed 
fiaring magnetic loop. 



5.1. SED Categorization Criteria 

As discussed above, in general we found that i?trunc correlates with Mdisk in the model SEDs, e.g., a 
small i?trunc can fit even a bare photosphere SED if M^^sk is made sufficiently small (see Q ■ Thus we found it 
helpful to visualize i?trunc versus Mjisk, as shown in the lower panel of Figs.[3]-in]in order to better interpret 
the SED model fits. A vertical line at 10"'^ M© indicates our disk mass threshold (see §1]). In orange, we 
show the stellar photosphere for reference. In each figure, the Mdisk and i?disk values corresponding to each 
model SED from the upper panel are shown as diamonds. For comparison, the magnetic loop height (i?ioop) 
and its uncertainty (see §2.ip are shown as a dashed line and a hatched region, respectively. 

From these plots, we can thus begin to assess the degree of spatial correspondence between i?trunc and 
-Rioop for each star. Cases for which i?trunc < -Rioop can be interpreted as representing disks that are within 
reach of the observed flaring loops. Cases for which i?trunc > ^loop are somewhat less straightforward to 
interpret because we must first account for the effects of dust sublimation. The location of dust destruc- 
tion, or the sublimatio n radius, is calculated using each model's stellar temperature and radius as follows 



( Robitaille et al.ll200d ): 



/1600K\ ^-^ 

-Rdust = -Rstar X I -7^: I (1) 

\ J star / 

where 1600 K is the dust sublimation temperature. Uncertainties in i?star and Tgtar create a range of possible 
values for i?dust- We adopt 5% uncertainty in Tstar and 20% uncertainty in -Rstar in this calculation (see 
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e.g.. lHillenbrandlll997l ). In each of Figs.[3H6](and in figures 3.1-3.32, accessible online), filled (red) diamonds 
represent models effectively truncated at the dust destruction radius, i?dust, while unfilled points have vertical 
bars in the —y direction to show where that particular model's dust destruction radius is located. For the 
model uncertainty in i?trunc, we adopt 5% and 50% for stars with Tcff <5200 K and Tcff >5200 K, respectively 
(see ^. Thus, when we say the disk is truncated "effectively" at the dust destruction radius, the intended 
meaning is that i?trunc and i?dust are equivalent within their uncertainties. We also report fiducial i?dust 
values for our sample stars in Table [TJ these are calculated using Eq. [T] and the i?star and Tstar data in that 
same table. 

The broadband fluxes used here trace the spatial extent of a disk's dust; in principle the inner edge of 
the gas in the disk could extend even closer to the stellar surface. For the cases in which we find that an 
observed dust disk is truncated at the dust destruction radius, the dust disk is likely truncated by sublimation, 
a process which would not remove gas. Indeed, some systems hav e been observed to be accreting even though 
the dust disk is truncated far from the star (JEisner et al.ll2005l l2007l) . For dust disks truncated near the 
sublimation radius, but not within reach of the magnetic loop, it is possible that a gas disk extends closer to 
the s tar and is truncated within the loop height (this has been observed bv lNaiita et al.ll2003t Eisner et al. 



20051 ). Conversely, if the dust disk is truncated outside the dust destruction radius, some other process may 
be responsible for cleari ng out the inner p ortion of the disk, and therefore we assume that the inner gas is 
cleared out as well (e.g.. llsella et al.ll2009l ). 



Finally, for a few stars we lack sufficient photometric data to adequately constrain the location of i?trunc. 
In most cases, this is due to a lack of Spitzer photometry and thus the longest wavelength measurement is 
the 2mass 2.2 /zm flux. Consequently, the model SEDs in these cases are largely unconstrained and result in 
a wide variety of possible star-disk configurations which can fit the observed SED. 

Based on these considerations, in what follows we categorize our sample stars into four groups, based 
on the degree to which the SEDs indicate that the inner disk edge is within reach of the flare loop height: 

Category 1: i?tiunc < ^loop^ The inner disk edge is clearly within reach of the magnetic flaring loop. 

Category 2: i?trunc > -Rioop but i?trunc ~ -Rdust^ The dusty inner disk edge is beyond the flaring loop 
height, however the dust disk is truncated at the dust-destruction distance and thus a gas disk may 
extend inward to i?ioop (i.e., i?trunc < -Rdust)- 

Category 3: i?trunc > -Rioop and i?trunc > ^dust: The inner edge of the dust disk is clearly beyond reach of 
the magnetic flaring loop. 

Category 4: Indeterminate: More than one category above is permitted by the available data (generally 
due to lack of Spitzer data). 



Recall that we take half the loop length as the flare loop height, as a conservative upper limit (see ij2.ip . 
In general, for a given star there are multiple SED models that are good flts to the observed SED, and in 
some cases the multiple best-fitting model SEDs yield a mixed verdict regarding the placement of i?trunc with 
respect to i?ioop. Thus if one of categories 1-3 above is favored by more than | of the best-fit BED models, 
we assign the star to that category, and we assign "indeterminate" (category 4) otherwise. 
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5.2. Example Cases 

As an example of our approach to interpreting the SEDs of our study sample, we show in Fig. [3] (upper 
panel) the SED of COUP 1410. The fiducial SED model in Fig. |31 corresponding to a disk with -Rtiunc = ^loop, 
predicts an excess of IR flux at wavelengths as short as 3 /im, unlike the data and best-fit model SEDs which 
follow the profile of a bare stellar photosphere to 4.5 /xm. Intuitively, this implies that the best-fit model 
SEDS must therefore correspond to disks with moderately large inner holes. Indeed, the lower panel of Fig. [3] 
shows that nearly all of the best-fit model SEDs, representing disks with 10~^ < A/disk/M© < 10^^, have 
-Rtrunc > 1 AU. Furthermore, the majority of these models are truncated well outside their respective dust 
destruction radii (i.e., i?trunc > ^dust); only one fit model has i?trunc ^ ^dust, and this model has very low 
Mdisk, below our threshold of 10^"^ Mq 

Note that the observed SED for this star does not in fact require any disk at all; the fact that many of the 
SED models shown in Fig. [3] exhibit large excesses longward of 4.5 /im implies only that these hypothetical 
disks with very large inner holes are formally permitted by the available data. These models thus provide a 
lower limit to the size of i?trunc that any as- yet undetected disk could possibly have. Since this lower limit is 
in this case much larger than -Rioop: we conclude that no disk is present that could interact with the observed 
magnetic flaring loop, and we assign COUP 1410 to category 3 ([J57 



COUP 141 (Fig. |4|) is a case in which the observed SED is reasonably well matched by the fiducial SED 
model, for which i?trunc = ^loop- The best-fit SEDs have inner truncation radii well beyond reach of the 
magnetic loop. However, these models' inner disk radii are also equal to their dust destruction radii, and 
thus it is likely that sublimation is responsible for the apparent clearing of the inner disk. In cases like 
COUP 141, while the magnetic loop may not intersect the dust disk, it could nonetheless intersect a gas disk 
that extends inward of the dust to within reach of i?ioop- Indeed, both Ca II and A(C/ — V) strongly indicate 
active accretion (Table [ij. Thus we assign COUP 141 to category 2. 

As another example, consider COUP 720 (Fig. [S|) . In this case, all of the best- fit model SEDs with Aidisk 
above our adopted threshold of 10^'^ M© have similar i?trunc ~ 0.1 AU, which overlaps -Rioop within its 
uncertainty. Many of these best- fit models, moreover, have i?trunc ~ ^dust, and thus may possess gas disks 
that extend even closer to the star. COUP 720 thus represents a good example of an SED that is consistent 
with -Rtrunc ~ ^loop, and for which the large magnetic loops observed by COUP may facilitate the magnetic 
star-disk interaction envisaged in magnetospheric accretion models. COUP 720 is assigned to category 1. 

Finally, consider COUP 997 (Fig. [6l). The observed SED data (0.34-4.5 /im) show excess IR flux. About 
half of the best-fit models are truncated at their dust destruction radii (category 2), while the other half 
are truncated beyond 1 AU (category 3). We also do not have Ca II or A(C/ — V) measurements to help 
disambiguate the two possibilities, and thus it is not possible to say which set of models correctly describes 
the observed star-disk system. Requiring additional data (particularly longward of ^ 10/im) to discriminate 
between the category 2 and 3 model fits, we assign this object to category 4. 



6. RESULTS 

In Table [SJ we present a summary of the results for the 32 stars in our sample. Following the procedure 
described in fj5l we have identified which stars' SEDs have massive (i.e., >10~^Mo) disks that are consistent 
with being within reach of the observed magnetic loops. Notes about each star relevant to its classification 
are provided in the figure captions. 
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We find six stars that appear to have SEDs consistent with -Rtrunc < ^loop (category 1). Another four 
stars do not show direct evidence of disks within reach of the magnetic loops, but could potentially have 
gas that extends interior to the observed dusty inner edge of the disk (category 2). Fourteen stars either 
have disks whose inner edges are situated beyond the reach of the magnetic loops, or are simply devoid of 
detected disk material entirely (category 3). For eight stars, we could not assign a definitive category as 
additional data are necessary to support or eliminate difi'erent classes of best-fit SED models (category 4). 



7. DISCUSSION AND CONCLUSIONS 

Of the 24 stars in our 32-star sample for which we have enough optical-infrared data to constrain the 
location of the disk inner edge (i.e., excluding stars in category 4; i j5.ll and Table [6]), for about 58% we are 
able to rule out close-in, mas sive disks within rea ch of the observed large flaring loops (category 3). These 
energetic flares discovered by iFavata et al.l (|2005i ) are evidently intrinsically stellar phenomena. This gives 



added justiflcation a posteriori for the application of the solar-flare cooling loop model to these stars ( ij2.ip . 
and suggests that it may be possible by further extension of the solar analogy to infer other flare-related 
properties for these flares, such as coronal heating rates and coronal mass ejections. The latter in particular 
may be important for furthering our understanding of mass and angular momentum loss in these low-mass 
PMS stars. 

Our sample also includes six cases for which the SED clearly indicates a dusty disk that extends close 
enough to the star to permit interaction with the flaring loop (category 1; Table [6]). In four additional 
cases the dust disk appears to be truncated beyond the reach of the flaring loop, but at or close to the 
predicted dust destruction radius (category 2). In these cases, the dust disk may in fact be truncated by 
dust sublimation, a process which does not remove gas. Thus, in category 2 objects, it seems likely that 
a gas disk (undetected in the broadband flux measurements used in our SED models) extends closer to the 
star and may be within reach of the observed flaring loops. One of the four category 2 objects has a Ca II 
measurement from the literature (couP 141), and two have A{U — V) measurements (couP 141 and COUP 
1568). COUP 141 is, interestingly, the most strongly accreting object in the sample as probed by its Ca II 
equivalent width (see Table[l]) and it also has a very negative A(C/ — V); the combination of these indicators 
is strong evidence for ongoing accretion. For COUP 141 the AKs near-iR excesses also indicates a disk very 
close to the star (see Sec. [3] and Table |4]). These examples further strengthen the interpretation of the 
category 2 stars as likely having gas accretion disks within reach of the stellar magnetosphere. 

Several studies clearly find a high f reque ncy of close-in, massive dusty disks in the ONC population as a 



whole. For example, iHillenbrand et al.l (|1998l ) find a disk fraction in the ONC of ~70% on the basis of excess 
emission at 2.2 /im. Thus, the ~25% (category 1) or ~38% (categories 1 and 2) frequency of close-in disks 
in our 32-star sample, representing the ~1% of ONC stars with the most powerful X-ray flares observed by 
COUP, is evidently not representative of the disk characteristics of the ONC as a whole. 

The six cases in which we have observed stars with close-in, dusty disks which intersect the magnetic 
loops (category 1) are interesting candidates for further study. Specifically, it would be informative to 
determine if the flares in these cases are in some way different than the category 3 cases. Three of these 
objects have A{U — V) measurements, all of them < —0.3 (Table [T]), strongly indicat ing active accretion . 



Temporally linking accretion as seen in optical variability to the X-ray flare events (e.g.. lStassun et al.ll2006l ) 
could solidify whether a magnetic star-disk interaction has taken place. Geometric information would be 
necessary to determine where the magnetic loop is on the stellar surface — for example, a given flaring loop 
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could extend in a direction perpendicular to the disk and thus not interact, even if the dust (or gas) disk is 
within the appropriate distance from the stellar surface. 

The question of how the large magnetic s tructures are stabil ized (prior to the flaring event) was posed 
in the discovery publication of these objects (JFavata et alj 120051 ). It was proposed that the loops may be 
anchored to corotating disk material and thereby not subject to shear which could disrupt the loops. In 
this work we have found such disk- anchoring to indeed be a possibility for 10 objects (categories 1 and 2), 
but further analysis is needed to determine how massive and ionized a disk must be to enable interaction. 
For fourteen objects (category 3), disk-supported loops are unlikely, as these stars lack massive disks within 
reach of the loops. We speculate that as long as the confining magnetic field at the upper end of the loop 
is suffici ently strong and the confined material corotating, stability is feasible even without a disk. For 
example. ICranmeii (2009!) describes a loop geometry in which the pressure of the confined gas decreases with 
increasing loop length, implying that in fact the largest coronal loops may be most stable against rupture. 
Additionally, our findings may imply that the largest fiaring loops cannot readily form in the presence of a 
disk, given that they appear in our sample to preferentially occur on stars lacking close-in disk material. 

Alternate scenarios to expla in the long X-ray decay timescales inc lude observations of corotating, embed- 
ded structures in coronae (e.g.. Collier Cameron &: Robinsonlll989ai br). and prominences and stellar winds 
( Skellv et al.ll2008t iMassi et al.ll2008[) . The observed coronal structures from those studies lend additional 
support to the idea of these large magnetic structures remaining stable in hot coronae or within the stellar 
wind over multiple rotation periods as t he X-ray flare decays, even if no disk is present to anchor the magnetic 



loop (jjardine fc van BallegooijenlbOOSl ) 



In summary, the 32 most powerful fiares observed by the COUP survey were found to h ave magnetic 



struc tures multiple stellar radii in arc length confining the x-ray emitting, heated plasma (JFavata et al. 



20051 ). With the goal of understanding the nature of these large x-ray emitting flare structures, we have 
modeled the optical-infrared SEDs of these objects, finding 58% to be lacking close-in circumstellar disks to 
which these loops could anchor. It is evident that in at least these cases the large-scale flares are phenomena 
of purely stellar origin, neither triggered nor stabilized by star-disk interactions. 
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Table 1. Stellar Parameters 



Object 
Name 



Toft" 

fKl 



Mass" 
[Mq] 



Radius^ 

[R©] 



Magnetic Loop 
Height [Rq]" 



Dust Destruction 
Radius*^ [R©] 



Ca II EW 

[A] 



Literature 
Av" [mag] 



Re-calculated 
Av" [mag] 



[mag] [mag] 



COUP 7 


4581 


2.12 


6.23 




0.65 


56.8 




0.75 


0.67 


-0.58, 


+0.54) 


0.12 


0.03 


COUP 28 


3802 


0.53 


2.3 


7.91±5.46 


14.2 


1.6 


0.63 


0.30 


-0.30, 


+0.54) 


0.20 


0.24 


COUP 43 


3606 


0.4 


2.92 


16.1 


(+3.16,-5.32) 


16.0 


1.4 


1.36 


1.18 


-0.60, 


+0.56) 


0.012 


0.64 


COUP 90 


3802 


0.52 


2.51 


1.04 


(+4.86,-1.04) 


15.5 


1.6 


4.97 


3.97 


-1.13, 


+ 1.13) 


0.053 




COUP 141 


5236 


2.11 


3.3 




1.97 


39.2 


-17.8 


1.83 


2.17 


-0.63, 


+0.60) 


0.37 


-0 


67 


COUP 223 


4395 


1.19 


2.79 




6.95 


23.3 


1.7 


4.66 


5.81 


-0.75, 


+0.88) 


0.31 






COUP 262 


4395 


1.13 


1.58 


28.5 


(+16.8,-22.3) 


13.2 


2.3 


3.77 


7.89 


-1.20. 


+ 1.24) 


0.14 






COUP 332 


3111* 


0.5* 


2* 




L05+38.4 


8.09 






12.8 


-1.67, 


+2.14) 


0.032 






COUP 342 


4729* 


0.5* 


2* 


20.3 


(+21.4,-20.3) 


19.5 






7.89 


-0.78, 


+0.85) 


0.085 






COUP 454 


4775 


2.35 


4.58 


46.4 


(+14.8,-9.20) 


45.5 


2.1 


5.85 


6.39 


-0.85, 


+0.87) 


0.077 






COUP 597 


5662 


1.49 


2.01 


3.16 


(+7.62,-3.16) 


28.6 


4.5 


2.69 


3.32 


-1.35, 


+ 1.39) 


0.27 


-0 


26 


COUP 649 


3589 


0.4 


2.17 




9.20 


11.8 





4.11 


3.82 


-0.73, 


+0.83) 


0.17 






COUP 669 


4581 


1.52 


2.59 


13.2 


(+2.73,-3.31) 


23.6 




1.96 


2.33 


-1.10, 


+ 1.20) 


0.049 






COUP 720 


4452* 


0.5* 


2* 


19.0 


(+69.3,-19.0) 


17.2 






11.8 


-1.33, 


+ 1.66) 


-0.23 






COUP 752 


3802 


0.54 


1.67 


9.35 


(+1.29,-9.35) 


10.3 


1.1 


0.07 


0.64 


-0.64, 


+ 1.10) 


0.047 


-0 


70 


COUP 848 


3342 


0.29 


1.98 


23.3 


(+4.03,-3.31) 


9.31 





1.72 


1.35 


-1.34, 


+ 1.46) 


0.093 






COUP 891 


4775 


2.43 


4.85 


24.9 


(+3.88,-3.31) 


48.8 


1.8 


8.00 


10.7 


-1.02, 


+ 1.12) 


0.13 






COUP 915 


4613* 


0.5* 


2* 


11.2 


(+3.88,-2.59) 


18.4 






15.8 


-1.76, 


+2.56) 


-0.065 






COUP 960 


3177 


0.24 


2.16 


0.53 


(+1.48,-0.53) 


9.27 





2.72 


1.29 


-0.41, 


+0.43) 


-0.013 






COUP 971 


3999 


0.69 


3.28 


5.03 


(+0.72,-0.43) 


22.5 


1.8 


0. 


0.00 


-0.00, 


+0.15) 


-0.24 


-0 


80 


COUP 976 


3177 


0.18 


0.91 




10.9 


3.85 





0. 


2.88 


-0.69, 


+ 1.13) 


0.28 






COUP 997 


3856* 


0.5* 


2* 


4.89 


(+4.03,-3.88) 


12.7 






1.82 


-1.14, 


+ 1.13) 


0.060 






COUP 1040 


4281* 


0.5* 


2* 


1.25 


(+1.91,-1.25) 


15.8 






16.5 


-2.03, 


+3.53) 


-0.10 






COUP 1083 


4698* 


0.5* 


2* 


33.8 


(+9.63,-7.05) 


19.2 






4.16 


-0.56, 


+0.61) 


-0.0072 






COUP 1114 


4903* 


0.5* 


2* 


9.78 


(+2.30,-9.78) 


21.0 


-1.5 




6.08 


-0.71, 


+0.74) 


0.086 






COUP 1246 


3177 


0.23 


1.62 


5.75 


(+1.01,-1.15) 


6.95 





0.92 


1.52 


-0.60, 


+0.70) 


0.45 






COUP 1343 


3649* 


0.5* 


2* 


13.8 


(+2.45,-2.44) 


11.6 






3.04 


-0.83, 


+3.66) 


0.60 






COUP 1384 


3802 


0.52 


2.46 


7 


.33+4.46 


15.1 


1.9 


0. 


0.67 


-0.67, 


+0.82) 


-0.095 


-0 


35 


COUP 1410 


3606 


0.36 


0.51 


15.8 


(+12.9,-15.8) 


2.80 





0.57 


4.98 


-0.98, 


+0.97) 


0.053 






COUP 1443 


5528* 


0.5* 


2* 


3.74 


(+4.17,-0.57) 


13.6 






0.26 


-0.26, 


+0.55) 


0.18 


0.14 


COUP 1568 


5236 


2.55 


3.99 


0.53 


(+1.62,-0.53) 


48.1 




0.59 


1.06 


-0.56, 


+0.60) 


-0.040 


-0.22 


COUP 1608 


3724 


0.48 


1.76 


11.8 


(+2.44,-11.8) 


10.8 


-1.3 


0.93 


0.25 


-0.25, 


+0.41) 


0.43 


-1.2 



I 



^Taken from JHillenbrandl j 19971) , unless otherwise noted. 

Taken from lFavata et al.| 1 120051) . No uncertainty is quoted in cases where lFavata et al.| ||2005|) used only two points in fitting the flare log(T)-log{ne) decay slope. 
■^Derived in this work. 



Note. — t Temperatures not in literature; assigned T^ff from best fit SED model. | Fiducial masses and radii of 0.5 Mq and 2 Rq taken in cases for which there 
were no measurements available in the literature. 



Table 2. Hubble ACS and WFI Fluxes 



Object 
Name 



4317A 
Flux [injy] 



5359A 

Flux [inJy] 



6584A ^^^ 
Flux [mjy] 



7693A 
Flux [mJy] 



9055A 
Flux [mJy] 



0.36/iin 
Flux [mJy] 



0.44/j.m 
Flux [mJy] 



0.55/im 
Flux [mJy] 



0.83fj.ni 
Flux [nijy] 



COUP 7 






170 






3.9±0.36 


40. ±3. 7 


83.0±7.6 


290±27 


COUP 28 


0.89±0.082 


2.9±0.27 


7.3 






0.099±0.0091 


1.1±0.10 


2.7±0.25 


16±1.5 


COUP 43 


0.40±0.037 


1.6±0.15 


5.2 


10. ±0.95 


20. ±1.9 


0.033±0.0031 


0.42±0.039 


1.2±0.11 


15±1.4 


COUP 90 


0.028±0.0026 


0.15±0.014 


0.69 


1.7±0.15 


4.0±0.37 




0.038±0.0035 


0.13±0.012 


2.9±0.26 


COUP 141 


7.9±0.73 


16±1.4 








1.3±0.12 






66±6.1 


COUP 223 


0.074±0.0068 


0.32±0.029 


1.6 








0.095±0.0088 




6.2±0.57 


COUP 262 


O.OlliO.OOlO 


0.095±0.0088 


0.57 


2.0±0.18 


6.7±0.61 






0.11±0.010 


4.4±0.41 


COUP 332 








0.049±0.0045 


0.33±0.030 










COUP 342 


0.008±0. 00072 


0.078±0.0072 


0.52 


1.9±0.17 


6.3±0.58 






0.099±0.0091 


4.4±0.41 


COUP 454 


O.lliO.OlO 


0.64±0.059 


2.6 


6.8±0.62 


16±1.5 




0.16±0.014 


0.63±0.058 


13±1.2 


COUP 597 


2.0±0.18 


5.8±0.54 






26±2.4 


0.34±0.032 


2.3±0.22 


5.3±0.48 


23±2.1 


COUP 649 


0.019±0.0017 


0.10±0.0093 


0.38 


1.5±0.14 


3.4±0.31 








2.7±0.25 


COUP 669 


1.8±0.17 


6.6±0.61 


15 








2.4±0.22 


5.6±0.52 


31±2.8 


COUP 720 






0.070 


0.31±0.029 


1.5±0.14 








0.99±0.091 


COUP 752 


0.85±0.078 


2.6±0.24 


7.1 


10. ±0.92 


16±1.4 


0.18±0.016 


1.3±0.12 


3.0±0.27 


16±1.5 


COUP 848 


0.097±0.0089 


0.37±0.034 


1.5 


3.3±0.30 


6.9±0.63 






0.46±0.042 


5.4±0.50 


COUP 891 


0.002±0. 00022 


0.046±0.0043 


0.40 


2.1±0.20 


9.1±0.84 






0.046±0.0042 


5.7±0.53 


COUP 915 


















0.17±0.015 


COUP 960 


0.022±0.0020 


0.082±0.0075 


0.40 


1.6±0.14 


4.3±0.40 




0.029±0.0027 


0.076±0.0070 


3.1±0.28 


COUP 971 


10. ±0.96 


21±2.0 


43 






1.6±0.15 


11±1.0 




64±5.9 


COUP 976 


0.016±0.0015 


0.080±0.0074 


0.51 


1.4±0.13 


3.8±0.35 








2.7±0.25 


COUP 997 


0.60±0.055 


2.1±0.19 


5.9 


11±0.98 


18±1.7 




0.77±0.071 


2.0±0.18 


14±1.3 


COUP 1040 


















0.10±0.0093 


COUP 1083 


0.43±0.040 


1.8±0.17 


6.2 


14±1.3 


26±2.4 




0.55±0.051 


1.7±0.16 


18±1.7 


COUP 1114 


0.23±0.021 


1.6±0.15 


7.0 








0.38±0.035 


1.6±0.14 


37±3.4 


COUP 1246 


0.046±0.0043 


0.22±0.020 


1.2 


1.3±0.12 


5.2±0.48 




0.059±0.0055 


0.19±0.017 


3.9±0.36 


COUP 1343 


0.058±0.0054 


0.31±0.028 


1.7 


2.8±0.26 


6.7±0.62 




0.082±0.0075 


0.27±0.025 


5.0±0.46 


COUP 1384 


2.8±0.26 


7.9±0.73 


17 






0.39±0.036 


3.3±0.30 


6.9±0.64 


30. ±2. 7 


COUP 1410 


0.0040±0. 00037 


0.022±0.0020 


0.13 


0.48±0.045 


1.4±0.13 






0.021±0.0019 


1.0±0.093 


COUP 1440 


1.9±0.18 


6.2±0.57 


15 




24±2.2 


0.21±0.019 


2.5±0.23 


6.0±0.56 


31±2.8 


COUP 1568 




85±7.8 








8.4±0.77 


50. ±4. 6 




202±19 


COUP 1608 




1.4±0.13 








0.20±0.018 


0.73±0.068 




11±1.02 


Zero points 


25.793'^' 


25.744<^' 


22.393'^' 


25.291'^' 


24.347'^* 


1823 Jy'^' 


4130 Jy<-'' 


3640 Jy*^' 


2430 Jy'^' 



I 



Note. — The first five data eolumns report HST fluxes, the last four columns report WFI fluxes. Zcropoint fluxes used are reported in the final row of the table. 

^ 'For the purposes of fitting, the Ha fluxes were given a wide (99%) error in order to allow for variability commonly observed T Tauri stars. 

^ 'The ACS data utilized were in the Vegamag system; these values were used in the conversion to the ABniag system, for which fiux calculation is straightforward 
as outlined in the online ACS documentation found here: |http: //www. stsci.edu/hst/acs/aiialysis/zeropoints/| 

^ 'Johnson-Cousins zero points; UBVIc- 



Table 3. V, /, 2MASS, and Spitzer Fluxes 



Object 
Name 



0.55/^iin 
Flux [mjy] 



0.79;/m 
Flux [mJy] 



1.235 /j-in 
Flux [mJy] 



1.662 /iin 
Flux [mJy] 



2.159 Mm 

Flux [mJy] 



3.6 ^m 
Flux [mJy] 



4.5 /^in 
Flux [injy] 



5.8 /j.ni 
Flux [nijy] 



8.0 ^m 23.6 ^m 

Flux [mJy] Flux^ [mJy] 



COUP 7 


102±9.4 


280±25 


460±42 


590±54 


440±41 


200±18 


130±12 


84±7.8 


50. ±4. 6 


10 


COUP 28 


3.8±0.35 


17±1.6 


39±3.6 


47±4.4 


41±3.8 


20. ±1.8 


13±1.2 


8.9±0.82 


6.5±0.60 


14 


COUP 43 


2.2±0.20 


15±1.4 


51±4.7 


72±6.6 


62±5.7 


33±3.0 


23±2.1 


16±1.4 


11±0.98 




COUP 90 


0.084±0.0078 


1.8±0.16 


14±1.2 


27±2.5 


27±2.5 


16±1.5 


11±1.0 


6.6±0.60 






COUP 141 


24±2.2 


71±6.6 


130±12 


180±16 


170±16 


150±14 


140±13 


120±11 


160±15 




COUP 223 


0.42±0.038 


5.1±0.47 


39±3.6 


93±8.6 


120±11 


140±13 


130±12 


110±10. 


91±8.4 




COUP 262 


0.31±0.028 


2.7±0.25 


35±3.2 


96±8.8 


130±12 


98±9.0 


68±6.3 


39±3.6 






COUP 332 






4.8±0.44 


22±2.0 


40. ±3. 7 


64±5.9 


69±6.3 


87±8.0 


150±14 




COUP 342 






32±3.0 


81±7.5 


95±8.7 


72±6.6 


50. ±4. 6 


63±5.8 






COUP 454 


0.66±0.061 


10. ±0.93 


74±6.8 


150±14 


150±14 


150±14 


130±12 


120±11 


150±14 




COUP 597 


6.1±0.56 


21±1.9 


41±3.8 


58±5.4 


63±5.8 


71±6.5 


62±5.7 








COUP 649 


0.082±0.0075 


1.6±0.15 


14±1.3 


29±2.6 


31±2.8 












COUP 669 


5.8±0.54 


24±2.2 


68±6.3 


97±8.9 


83±7.7 


110±9.9 










COUP 720 






13±1.2 


43±4.0 


65±6.0 


83±7.7 


70±6.5 


58±5.4 


57±5.2 




COUP 752 


3.4±0.31 


12±1.1 


31±2.8 


40±3.7 


34±3.1 


20. ±1.9 


21±1.9 


15±1.4 


20. ±1.8 


11 


COUP 848 


0.36±0.033 


4.1±0.37 


17±1.6 


22±2.0 


20. ±1.9 












COUP 891 


0.10±0.0094 


3.4±0.31 


69±6.3 


240±22 


350±32 


240±22 


170±15 


120±11 


79±7.3 


90 


COUP 915 






5.6±0.52 


34±3.1 


68±6.3 


78±7.2 


67±6.2 


56±5.1 






COUP 960 


0.093±0.0086 


2.1±0.19 


llil.O 


13±1.2 


11±1.0 


6.6±0.61 


4.7±0.43 






36 


COUP 971 


26±2.4 


62±5.7 


102±9.4 


140±13 


82±7.6 


44±4.1 










COUP 976 






13±1.2 


24±2.2 


26±2.4 


24±2.2 


20. ±1.8 








COUP 997 


2.0±0.19 


11±0.98 


36±3.3 


55±5.0 


50.±4.6 


39±3.6 


27±2.5 








COUP 1040 






4.0±0.37 


26±2.4 


58±5.3 


75±6.9 


76±7.0 


60. ±5. 6 


56±5.2 




COUP 1083 


3.2±0.29 


16±1.5 


57±5.2 


95±8.8 


86±7.9 


54±5.0 


36±3.3 








COUP 1114 


1.9±0.17 


26±2.4 


190±18 


380±35 


410±37 


250±23 


190±18 


130±12 


86±7.9 




COUP 1246 


0.28±0.025 


3.2±0.29 


15±1.3 


24±2.2 


28±2.5 


18±1.7 


17±1.5 


21±1.9 


42±3.9 




COUP 1343 


0.32±0.029 


3.5±0.32 


25±2.3 


58±5.4 


82±7.5 


70±6.5 


62±5.7 


35±3.2 






COUP 1384 


8.2±0.75 


28±2.6 


66±6.1 


87±8.0 


69±6.4 


34±3.1 


23±2.1 


18±1.6 






COUP 1410 


0.14±0.013 


0.72±0.067 


5.8±0.53 


12±1.1 


12±1.1 


7.2±0.67 


4.8±0.44 








COUP 1443 


6.5±0.60 


25±2.3 


57±5.3 


77±7.1 


62±5.7 


32±3.0 


20. ±1.9 


18±1.7 






COUP 1568 


109±10. 


207±19 


290±26 


300±27 


240±22 


130±12 


100. ±9. 6 


97±8.9 


130±11 


156±14 


COUP 1608 






26±2.4 


39±3.6 


46±4.2 


45±4.1 


46±4.3 


34±3.2 


43±3.9 


119±16 


Zero points [Jy] 


3640*^' 


2490'^* 


1594<3' 


1024*^' 


666.7'^' 


280. 9^''* 


179.7(^) 


115.W 


64. IS*'*' 


7.17(5) 



05 
I 



^Fluxes reported without error arc 3cr upper limits. 
^iBessei] Ji97i). 
'' jCousinsI | |1976| V 
^ ICohen et~ai] 120031) . 

Spitzer Science Center data calibration manual, accessible online at: |http : //ssc . spitzer . caltech . edu/ document s/cookbook/htttil/cookbook-node208 .html| 
^ jEngelbracht et al.l i2007t) . 
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Table 4: Reliability of Near-IR Excess as Tracer of Inner Disk Edge 





Rtiunc ^ Rdust 


J^tiunc ^ J^dust 


A{Ks) > 0.3 


141^ , 223, leOS'' 


1246*^ 


A{Ks) < 0.3 


332, 454, 720, 752, 976, 
1040, 1384, 1568 


7, 28, 43, 90, 262, 597, 
669, 891, 960, 971, 1083, 
1114, 1410, 1443 



"Object possesses Ca II in emission (Table [TJ a spectroscopic indicator of active accretion. 

''We consider the large AKs for COUP 1246 to be spurious (see discussion of this specific case in Q, so this object actually 
belongs in the lower right quadrant of the Table. 

Note. — Category 4 objects have been omitted, as a clear determination of -Rtrunc could not be made; see i|5.1l 



Table 5: Reliability of Near-UV Excess as Tracer of Inner Disk Edge 





Rtrunc ^ Rdust 


Rtrunc > Rdust 


A{U-V) <-0.3 


141^ , 752, 1384, 1608^ 


971 


A{U-V) >-0.3 


1568 


7, 28, 43, 597, 1443 



^Object possesses Ca II in emission (Table [TJ a spectroscopic indicator of active accretion. 

Note. — Category 4 objects have been omitted, as a clear determination of Rtrunc could not be made; see ^5.11 

Table 6. Spectral Energy Distribution Result Summary 



Category 1 


Category 2 


Category 3 


Category 4 


332 


141 


7 


223 


454 


597 


28 


342 


720 


1040 


43 


649 


752 


1568 


90 


848 


1384 






262 


915 


1608 






669 


976 










891 


997 










960 


1343 










971 












1083 












1114 












1246 












1410 












1443 





Note. — Categories as described in ^5. 11 for all COUP sam- 
ple objects. 
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Fig. 1. — SED for COUP 26 2. Diamonds are photom etric data from Tables [IHSl (see i)2.2l) . The purple square 
is the iT-band flux from ( Hillenbrand et al.l 119981) used to calculate AXs for the red atmosphere model. 
The green curve is the best fitting NextGen stellar atmosphere model with Teff set to the literature value 
(Tabled]). The best- fit extinction, Ay, is reported upper righ t. The red curve re presents the same stellar 
atmosphere model but with Ay as previously determined by iHillenbrandl (|l997n based on a fit to the V 
an d I fluxes on l y. Ve rtical lines indicate the wavelengths of the (from left to right) V and Ic bands used 
by IHillenbrandl (|l997l ) and the WFI /-band newly reported here. The resulting /S.Ka, color excesses for both 
model atmosphere fits are reported at upper right (see also Table |T]) . It is clear that the new fit to the full 
set of available photometric fluxes results in a more accurate representation of the stellar SED. Whereas this 
star was previously identified as possessing a very large near-iR excess, the new SED fit here indicates no 
significant excess. 
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Fig. 2. — Effects of changing surface gravity (logg) on flux in the Ks passband (represented by vertical 
dashed Unes) as a function of steUar temperature. Each pane represents a different temperature within the 
expected range for young, low-mass stars in our study sample. In each plot, solar-metallicity NextGen 
stellar atmospheres are plotted with six different log g values. Each atmosphere is normalized to the /wfi 
bandpass at 0.83 /xm, indicated by red diamonds. For cooler stars, i.e., the atmospheres in the upper panels, 
Ks flux varies by a factor of three depending on log g; the effect is most pronounced at T^s ^ 4000 K. 
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Fig. 3. — SED of COUP 1410, category 3. All models are consistent with category 3; the SED is consistent 
with a bare photospher e. This object is d iscuss ed in greater detail in ij5.2l Upper panel: Best fit model SEDs 
from the model grid of iRobitaille et al.l (|2007l black dash-dotted curves). The fiducial ttsre model with 
^trunc = ^loop (see g]), showu as a solid blue line and normalized to the peak near-iR flux (J, H, or Ks 
band), is meant to illustrate approximately how the SED would appear if there were a disk within reach of 
the flaring magnetic loop. The dashed curve is a solar-metallicity NextGen atmosphere model representing 
the stellar photosphere. Red diamonds are measured fluxes (Tables [2H31) as detailed in t j2.2l Lower panel: 
Comparison of i?ioop and i?trunc for best-fit SED models. A dashed line illustrates the loop height above 
the stellar surface (solid, orange). Uncertainty in the loop height is shown as a gray, hatched region. Open 
diamonds represent the -Rtrunc values of the best-fit SED models from the upper panel. For each model, a 
vertical bar indicates the location of -Rdust for that model (according to Eq. [T]). Filled red diamonds indicate 
models which have i?trunc ~ -Rdust (see ^for more detail). Finally, the vertical dash-dotted line indicates 
our disk mass threshold value, 10 ~^ M©; less massive disks do not represent the disks typical of T Tauri 
stars (see 35). Figures for all 32 sample stars are available in color electronically (Figs. 3.1-3.32). 
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Fig. 4. — SED of COUP 141, category 2. All symbols are as in Fig. [3] 2/3 models are consistent with category 
2. Excess flux in the IRAC bands indicates a dusty disk. Fourteen best-fit SEDs are plotted, but three of these 
are degenerate in Mdisk or i?trunc, representing four inclinations of the same star-disk configuration. Three 
sets of model SEDs have Afdisk > 10"'^ M©, and more than 2/3 of these are truncated at their respective dust 
destruction radii (red points in lower panel) . Thus these disks may possess gas that extends inward of i?dust 
to i?ioop- This object is very likely accreting, based on its strong Ca II emission (— 17.8A) and A{U — V) 
excess of —1.26. For further discussion of this object, see 
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Fig. 5. — SED of COUP 720, category 1. All symbols are as in Fig. [S] All 175 models are consistent with 
category 1. This object is discussed in detail in ^H 
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Fig. 6. — SED of COUP 997, category 4. All symbols are as in Fig. [3l Most of the best fit model results are 
split between categories 2 (11/20) and 3 (9/20), placing this object into category 4. This object is discussed 
in detail in i15.2l 
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